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The non-Ohmic effect of a high electric field on the out-of-plane magneto-conductivity of a lay- 
ered superconductor near the superconducting transition is studied in the frame of the Langevin 
approach to the time-dependent Ginzburg-Landau equation. The transverse fluctuation conductiv- 
ity is computed in the self-consistent Hartree approximation for an arbitrarily strong electric field 
and a magnetic field perpendicular to the layers. Our results indicate that high electric fields can be 
effectively used to suppress the out-of-plane fluctuation conductivity in high-temperature supercon- 
ductors and a signiflcant broadening of the transition induced by a strong electric field is predicted. 
Extensions of the results are provided for the case when the electric field is applied at an arbitrary 
angle with respect to the layers, as well as for the three-dimensional anisotropic regime of a strong 
interlayer coupling. 



I. INTRODUCTION 

Outside the critical region above Tc, in the absence 
of magnetic field and for small electric fields, the excess 
conductivity due to fluctuations of the superconducting 
order parameter can be explained by the Aslamazov- 
Larkinl theory, subsequently extended by Lawrence and 
Doniach^ for two-dimensional layered superconductors, a 
situation very much resembling the crystal structure in 
high-temperature superconductors (HTSC). In the pres- 
ence of a magnetic field, the fluctuation transport proper- 
ties of superconductors were initially treated^i^i^ in the 
non-interacting (Gaussian) fluctuation approach, which 
predicted a divergence at Tc{H) that is, however, not 
observed. The physical reason is the motion of vortices 
providing dissipation and hence a finite fiux-flow con- 
ductivity. Ikeda et ali and Ullah and Dorseji^ showed 
that the theoretical divergence can be eliminated by using 
the Hartree approximation, which treats self-consistently 
the quartic term in the Ginzburg-Landau (GL) free- 
energy expansion. This approach was applied for the 
longitudinal^-^ and Hall conductivity,^'^ as well as for 
the out-of-plane conductivity^ in the linear-response ap- 
proximation for a layered superconductor under perpen- 
dicular magnetic field. 

The fluctuation transport properties can be calculated 
in the linear-response approximation only for sufficiently 
weak electric fields that do not perturb the fiuctuation 
spectrumiii At reasonably high values of the electric 
field, the accelaration of the paired electrons is so large, 
that on a distance of the order of the coherence length 
they change their energy by a value corresponding to 
the fiuctuation Cooper pair binding energy»12. This re- 
sults in an additional, field dependent, decay mecha- 
nism, and leads to deviation of the current-voltage char- 
acteristics from the Ohm's law. In connection with the 
low-temperature superconductors, the non-Ohmic fiuc- 
tuation conductivity in the absence of magnetic field has 
been studied theoretically for the isotropic case^^'^^ and 
found experimentally on thin aluminum films^^ii^ For a 
layered superconductor the issue has been more recently 



addressed for the in-plane conductivity, starting from a 
microscopic approach^^ and subsequently in the frame of 
the time-dependent Ginzburg-Landau (TDGL) theory, in 
the Gaussiani^ as well as in the self-consistent Hartree 
approximation^ Several experimental investigations of 
the fluctuation suppression effect of high electric fields in 
HTSC were performed for the in-plane paraconductivity 
in zero magnetic field r^iSLSSiSiSi and a good agreement 
with the theoretical modelsiii^ was proven. 

The non-linear effect of a strong electric field under 
the simultaneous application of a perpendicular magnetic 
field on the in-plane fiuctuation conductivity and Hall ef- 
fect was recently addressed by the authors of the present 
pa,peir^^ for a layered superconductor in the Hartree 
approximation of the TDGL theory. It has been revealed 
that the simultaneous application of the two fields re- 
sults in a slightly stronger suppression of the supercon- 
ducting fluctuation conductivity, compared to the case 
when the flelds are applied individually, while the rela- 
tive suppression of the excess Hall conductivity turns out 
to be stronger than for the longitudinal one. Experimen- 
tal investigations of the fiuctuation suppression effect of 
strong electric fields under simultaneous application of a 
magnetic field are however lacking so far. 

Also the out-of-plane conductance of a layered super- 
conductor in the non-Ohmic regime of high electric flelds 
has been, to our knowledge, neither experimentally nor 
theoretically investigated up to the present. The pur- 
pose of this paper is to provide a theoretical approach 
to this issue, in the frame of the TDGL equation solved 
in the Hartree approximation, when both the magnetic 
and the electric field are applied perpendicular to the lay- 
ers. We shall thus be able to find the expression of the 
Aslamazov-Larkin (AL) contribution to the out-of-plane 
fiuctuation conductivity at any values of the magnetic 
and electric flelds, and predict for the AL term a signifl- 
cant supplementary suppression induced by a strong elec- 
tric field in layered HTSC. Based on the proportionality 
between the Cooper pair concentration and the density- 
of-states (DOS) part of the out-of-plane conductivity, 
known from the microscopic theory in the linear response 
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approximation,- • we shall also give an estimate of the 
non-Ohmic effect on the DOS contribution. The work 
will be completed by extending the calculation of the 
non-Ohmic fluctuation conductivity to the case when the 
electric field is applied at an arbitrary angle with respect 
to the layers. This calculation, particularly useful for in- 
vestigations on vicinal thin films, is necessary because 
the non-Ohmic current produced by a tilted electric field 
cannot be calculated as the superposition of the non- 
Ohmic currents produced separately by the in-plane and 
out-of-plane field components, as it was the case in the 
linear response approximation. Eventually, the results of 
the paper will be extended also to the three-dimensional 
(3D) anisotropic regime of a strong interlayer coupling. 



II. NON-OHMIC OUT-OF-PLANE 
FLUCTUATION CONDUCTIVITY 

In order to calculate the non-Ohmic out-of-plane con- 
ductivity under the presence of a magnetic field per- 
pendicular to the layers, we shall adopt the Langevin 
approach to the gauge-invariant relaxational TDGL 
equatioD^ii^ governing the critical dynamics of the super- 
conducting order parameter in the l-th superconducting 
plane: 
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where m and nic are effective Cooper pair masses in the 
a6-plane and along the c-axis, s is the distance between 
superconducting planes, and the pair electric charge is 
eo = — 2e. The order parameter ipi has the same physical 
dimension as in the three-dimensional case, and SI units 
are used. The GL potential a — a^e is parameterized by 
Y?l2mil = ^V2mc^gc and e = \n{T/To), with 



ao 



To being the mean-field transition temperature, while 
^0 and ^oc are, respectively, the in-plane and out-of- 
plane coherence lengths extrapolated to T = 0. The 
order parameter relaxation time is given by22i2S 

= nh^ / Idm^QkBT . The magnetic field B, perpen- 
dicular to the layers, is generated by the vector poten- 
tial in the Landau gauge A = {0,xB,0), with x and y 
the in- plane coordinates. Since we are interested in the 
out-of-plane conductivity, we consider the electric field E 
as being applied along the z-axis, and generated by the 
scalar potential ipi — —Esl. The Langevin white-noise 
forces Ci {x, y, t) that describe the thermodynamical fiuc- 
tuations must satisfy the fluctuation-dissipation theorem, 
(0 {x,y,t)Q, {x',y',t')) = 2T^^kBT5{x^x')6{y-y')5{t- 
t')Su'/s. 

The quartic term in the thermodynamical potential 
will be treated in the Hartree approximation, which 



results in a linear problem with a modified (renormalized) 
reduced temperature 



£ + &(|^,|')/ao. 



(2) 



The out-of-plane fluctuation conductivity Aazz will be 
eventually found by calculating the Josephson current 
density between the ^-th and I + l-th layers, which in the 
chosen gauge writes: 



and further 
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[(v^rvi+i) - (i^i^t+i)] , (3) 
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It is worth mentioning, at this point, that a slight mod- 
iflcation of the GL free-energy functional for the layered 
superconductors was recently proposed;^ which, besides 
the BCS-type Josephson couphng, allows for additional 
interlayer interactions that can contribute to the con- 
densation energy and give rise to energy savings that 
enhance Tq. The fluctuation spectrum of the proposed 
functional and consequently the fluctuation-induced ob- 
servables like in-plane paraconductivity and magneto- 
conductivity were found to be the same as for the usual 
Lawrence-Doniach free-energy, once the transition tem- 
perature of each bare layer is renormalized to its value 
enhanced through the interlayer energy savings. '^"'^ It can 
be easily verified that this fact is valid also for the out-of- 
plane fluctuation conductivity, since the additional terms 
in the GL functional of Ref. |^ do not involve the phase 
of the superconducting parameter, and consequently do 
not influence the deflnition Q of the transversal current 
density. The results of the present paper will be thus ap- 
plicable also for the kind of interlayer coupling proposed 
in Ref. 31. 

We proceed further by introducing the Fourier trans- 
form with respect to the in-plane coordinate y, the layer 
index I, and also the Landau level (LL) representation 
with respect to the x-dependence, through the relation: 



^iix,y,t) 



dk r'' dq^ ^ , , , 
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where the functions u„ (x) with n e N build the or- 
thonormal eigenfunction system of the harmonic oscil- 
lator hamiltonian, so that (— -I- Ae^B^x^) Un {x) = 
2heB {2n + 1) u„ {x). The TDGL equation (jTJ wiU write 
in the new variables: 



dt 



d_ 
dq 
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where the new noise terms (nik, q, t) are delta-correlated 
as {Uk,q,t)Q'{k',q',t')) = 2r-'kBT{2nf ■ 6{k - 
k')S{q - q')5{t - t')S„ 
notations: 



£„ = £ + (2n + 
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We have also introduced the 
B _ 2eilB 



r = 



26. 



h = 
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with h denoting the reduced magnetic field and r the 
anisotropy parameter. 

Equation © can be solved with the aid of the Green 
function technique and has the solution: 
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so that the correlation function between the order pa- 
rameter in two layers I and V will be given by: 
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where the electric field enters the parameter 
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with Eq = W^/SksT / Tre^o being the characteristic elec- 
tric field defined as in Refs. [l3and[li 

We point out that the sum over the LL in Eq. 
must be cut off at some index Nc, reflecting the 
inherent UV divergence of the Ginzburg-Landau the- 
ory. The classicali2i2£ procedure is to suppress the 



short wavelength fluctuating modes through a momen- 
tum (or, equivalently, kinetic energy) cut-off condi- 
tion, which, in terms of the LL representation writes2i2£ 
{2ehB/m) (?^ + ^) < cqq — ct? l2mQ^, with the cut-off 
parameter c of the order of unity. A total energy cut- 
off was also recently proposed,^- whose physical meaning 
was shown to follow from the uncertainty principle, and 
whose importance is revealed especially at high reduced- 
temperatures and magnetic fields close to 5^2(0)1^ For- 
mally, the total energy cut-off can be obtained from the 
momentum cut-off by replacing c with c—e. However, for 
low magnetic fields with respect to -Bc2(0) and in the crit- 
ical fluctuation region, the two cut-off conditions almost 
coincide quantitatively, so that we shall apply for sim- 
plicity the momentum cut-off procedure. In terms of the 
reduced magnetic field h, it writes thus h [Nc + |) = c/2. 

Now we are able to apply the expression ^ in Eqs. (0) 
and ^ in order to write, respectively, the self-consistent 
Hartree equation and the out-of-plane current density. 
Thus, after performing the q integral in the correlation 
function taken for I = I', one obtains the renormaliz- 
ing equation for the reduced temperature e, 



A', 



e = \n — + 2gThy" / e-"(^"+5)/, 



( — sin(pw) j , 
(11) 

where Iq{x) is the modified Bessel function and the 
parameter 



(12) 



was introduced according to the expression of the quartic 
term coefRcienlj^ b = iioKQ^CQfi^ /2m'^ , with kql being 
the in-plane Ginzburg-Landau parameter kql ~ Aq /Co- 
in an analogous manner, after computing the correla- 
tion function for I' = and using the current den- 
sity definition (|3J|, one can eventually obtain the out-of- 
plane fluctuation conductivity under arbitrary magnetic 
and electric fields, 



with /i(a;) the modified Bessel function of first order. The integrals in Eqs. Hll() and H13|l are convergent provided 
e -\- ^ -\- h > 0, so that e„ -I- | > for any LL index n. This condition is however assured while solving Eq. ()ll|l for 
the parameter e at any temperature T. 

It is useful to write Eq. H13|l also when the cut-off is neglected, i.e. for c, Nc — > 00: 



A.-— (i.,i.) = ^ ^^^f: . e-"(^-^i).g^/,(fsin(p.) 

Jo 1 — exp (— 2/iu) pu \2p 



(14) 
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This slightly simpler formula provides a good approx- 
imation for Eq. in the temperature region close to 
the transition (where e <C c) and for small magnetic fields 
(for which Nc — {c — h)/2h is already high). The cut-ofF 
procedure remains however essential for calculating the 

Cooper pairs density ^|'(/;;(x, t)|^^ contained in the renor- 

malization equation since the m- integral would be 

divergent for c, Nc ^ oo. As we shall see in Section IVlII 
Eq. I|14() can be also directly transformed in order to find 
the three dimensional limit (i.e. for s ^ 0) of the non- 

Ohmic fluctuation conductivity A(7^j"|^^^^^ parallel to 
the magnetic field, when the cut-ofF is neglected. 

In Eq. we have explicitly specified that the out- 
of-plane fluctuation conductivity Aa^^ corresponds to 
the Aslamazov-Larkin (AL) fluctuation process, since the 
phenomenological Ginzburg-Landau theory cannot ac- 
count for indirect contributions like the density-of-states 
(DOS) and Maki-Thompson (MT) terms, which can be 
found only from a microscopical approach. However, 
whereas the DOS and MT contributions to the in-plane 
paraconductivity and Hall effect coefficient are known to 
be negligible near the superconducting transition with 
respect to the AL term due to the more singular be- 
havior of the latter, an investigation of the out-of-plane 
conductivity needs taking into account also the DOS 
term, which can compete with the AL one especially for 
highly anisotropic materials, as pointed out by Larkin 
and Varlamovii^ The reason is that the DOS contribution 
to the out-of-plane fluctuation conductivity turns out to 
be proportional to a lower order of the interlayer trans- 
parency than the AL one, as shown by the microscopical 
approach in the linear response approximation&iS We 
shall tentatively give in Section IIVI an estimate of the 
DOS contribution to Aazz also for an arbitrarily strong 



electric field, after discussing the limit values of expres- 
sions and in the cases of a vanishing magnetic 
or electric field. 



III. LIMIT CASES B ^0 AND/OR E ^0 

The Ohmic out-of-plane fluctuation conductivity in the 
presence of a magnetic fleld but in an infinitesimally small 
electric field can be easily obtained by taking the limit 
p — > in Eq. H13I) , which acquires thus the form: 

2 2 ^<= 
n=0 

The expression (|15|) matches thus the result previously 
obtained within the diagrammatic microscopic approach, 
for Gaussian fluctuations (i.e. with 'e — e), in the 
linear response approximation^ Analogously, the self- 
consistent equation Hll|) becomes, in the same limit, 

T 

^E^o,B>o = ln— + 2gThJ2[enien + r)rKm 

which represents the Hartree renormalization equation 
in the linear response approximation under an applied 
magnetic field found in Ref. @. 

The other limit case, namely for vanishing magnetic 
field but under a finite applied electric field, needs taking 
the limit /i ^ in Eqs. and Hll|l after performing the 
sum over Landau levels and taking into account the cut- 
off condition, so that the Hartree self-consistent renor- 
malization and the out-of-plane fluctuation conduc- 
tivity l(T^ become, respectively. 
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(17) 
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It can be easily verified, by using the integral identities 
of the Bessel functions /q and Ii , that in the further limit 
E ^ 0, Eq. (O becomes 



~i 1 ^ I o Ve + c+Ve + c + r 
£]e=o,b=o = In ^^t + In , (19) 



as also found in Ref. \i% while expression (|18|l takes. 



respectively, the form: 



2^ l£;=o,s=o 



_v/F(i+7) 

£ + 0+ S 



V(e + c) (e + c + r) 



(20) 



If one neglects the cut-off procedure (i.e. c — > oo}, ex- 
pression H2U|) matches the result obtained in Ref. [13 for 
the AL out-of-plane fluctuation conductivity in the linear 
response limit and in the absence of magnetic field, with 
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the difference tliat in Ref. 0, based on tfie Gaussian ap- 
proximation, the reduced temperature e = In (T/Tq) is 
present instead of our Hartree renormahzed e. 

The Hartree renormalization procedure consists in us- 
ing the reduced temperature parameter e instead of e = 
In (T/To), by solving Eq. l(T^ . This procedure causes the 
critical temperature to shift downwards with respect to 
the bare mean- field transition temperature Tq. In anal- 
ogy with the Gaussian fluctuation case, we shall adopt as 
definition for the critical temperature the vanishing of the 
reduced temperature, e = 0, where the fluctuation con- 
ductivity, given by Eq. (|2U|I . diverges. In practice, one 
knows experimentally the actual critical temperature Tco 
measured at very low electrical field and with zero mag- 
netic field, so that the relationship between Tco and Tq 
will be found by putting e = in Eq. ifT^ll . It writes tiS 



cO 



(21) 



Now, having Tq one can use Eq. Hll|l for any temper- 
ature T and fields E and B in order to find the actual 
renormahzed £(T, B). 



IV. ESTIMATION OF THE NON-OHMIC DOS 
CONTRIBUTION 

In the microscopic approach of Ref. jQ, valid in the 
linear response approximation (i.e. for vanishing electric 
field) and for non-interacting (Gaussian) fluctuations, the 
DOS contribution to the out-of-plane fluctuation conduc- 
tivity under a magnetic field is found to amount, in our 
notations: 



in the — s- limit. The proportionality between the 
DOS fluctuation conductivity and the Cooper pair con- 
centration in Eq. is qualitatively easy to be grasped, 
since the DOS contribution means in fact the reduction 
of the normal state conductivity due to the decrease of 
the one-electron density of states, which reduction is in 
turn proportional to the superfluid density. 

We may assume that the proportionality (|24|) will hold 
also in the case of an arbitrarily strong electric field and 
in the Hartree approximation, so that the DOS contribu- 
tion to the out-of-plane fluctuation conductivity can be 
generally written: 
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where Eqs. ^ and are to be compared in order to 

reveal the Cooper pair density Analogous with 

Eq. H18|l we can consequently write the DOS contribution 
also in the vanishing magnetic field limit. 
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The equations (|26|l and (|27|l remain however to be 
confirmed or refuted by a microscopic approach. 



A DOS I Gaussian 



e S Kr 



TV, 



/i^[e„(en + r)] ^(^2) 



V. EXAMPLE: OPTIMALLY DOPED 

YBa2Cu306+^ 



where the parameter k depends on the impurity scatter- 
ing time r and temperature:^ 

(23) 

with ip (x) the Euler digamma function. One can notice 
that Eq. (|22|l can be written as: 



A DOS I Gaussian 



where 



2mckBT 



Gaussian 



E^O,B>0 



, (24) 



Gaussian 



£;^0,B>0 



mkBT 

TTh?S 



hJ2Wnien + r)] ^(25) 



is the Coooper pairs density for vanishing electric field 
and in the Gaussian approximation (i.e. with £ = e), as 
one can infer from the general correlation function 



In order to illustrate the main features of our model, 
we take as example a common HTSC material, like the 
optimally doped YBa2Cu306+x (YBCO). Typical char- 
acteristic parameters are then: s = 1.17 nm for the in- 
terlayer distance, = 1-2 nm and ^oc = 0.14 nm for 
the zero-temperature-extrapolated in-plane and out-of- 
plane coherence lengths, respectively, kql = 70 for the 
Ginzburg-Landau parameter, Tco — 92 K for the criti- 
cal temperature under very small electric and zero mag- 
netic field, and the parameter k = 3.57 that corresponds 
to a scattering time r ~ 30 fs in Eq. H23|) . It must 
be stated that the form of the normal-state background 
chosen for the temperature region masked by the onset 
of the superconductivity can be crucial for an eventual 
comparison with the experiment. There is however no 
consensus whether the peculiarities of the out-of-plane 
resistivity, namely its peak and its non-metallic charac- 
ter just above Tco, as observed in the oxygen-deficient 
YBCO and the more anisotropic Bi2Sr2CaCu208+2;, are 
mainly due to the competition between the fluctuation 
AL and DOS contributions, as illustrated in Ref. ^3 by 
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succesfull fits at different magnetic fields in the Ohmic 
regime wliile assuming a metallic linear extrapolation for 
the normal-state resistivity, or to the inherent behav- 
ior of the normal-state itself, as suggested by analysis 
of conductivity in incoherent layered crystals."^* Since in 
this paper we focus on the fluctuation conductivity in 
the non-Ohmic regime and need the normal-state back- 
ground only for illustration purposes of the resistivity 
characteristics, we shall further assume, for simplicity, 
an out-of-plane normal state resistivity almost constant 
near the transition, with a typical value = 4mi7cni 
for optimally doped YBCO^ 

In Fig. ^ the out-of-plane resistivity is presented, 
when both AL and DOS contributions are taken into ac- 
count, while Fig. ^ shows the effect of the AL term 
alone, if the same normal state background is assumed. 
One can notice the supplementary broadening of the 
transition induced by a strong electric field, and also the 
relative reduction of the non-Ohmic effect when a mag- 
netic field is simultaneously applied. The effect of various 
electric fields, at a fixed magnetic field, on the AL and 
DOS fluctuation conductivities is detailed in Figs. ^ and 
respectively. It turns out that the non-Ohmic effect 
is important only for the AL term, while the DOS one 
is little affected by an electric field of experimentally ac- 
cessible strength. This behavior stems from the peculiar 
dependence of Eq. (|26|l on the electric field only in the 
argument of the Besscl function. 

It is worth mentioning that for Bi2Sr2CaCu208+a;, for 
which the DOS contribution as such competes stronger 
with the AL one due to the higher anisotropy, the esti- 
mate of the non-Ohmic effect on the DOS term turns 
out however to be even more insignificant than for 
YBa2Cu306+a;, and thus almost undiscernable in the 
same range of electric fields. The reason is the much 
smaller anisotropy parameter r which reduces the effect 
of the Bessel function factor in Eq. (|2fi|l . 



arately by its components. However, this is not anymore 
the case in the non-Ohmic regime of a strong electric 
field. As we shall see below, the current components will 
depend now on all the field components, and not only to 
the particular one corresponding to the respective axis. 
This case requires therefore a special treatment, where 
both the in-plane and out-of-plane electric field compo- 
nents are to be included from the beginning in the TDGL 
equation. This calculation would be particularly useful 
if the investigation of the out-of-plane non-Ohmic con- 
duction were performed on vicinal thin films;^ where a 
mixture of the in-plane and out-of-plane transport prop- 
erties is assessed, since the injected current has a slanted 
direction with respect to the crystallographic axes. 

We shall consider in the following the case of an elec- 
tric field E applied on a layered superconductor at an 
angle 6 with the c-axis, having thus the components 
(^Ex — EsinO, Ey — 0^ Ez = Ecos9), and generated by 
the scalar potential tp — —E^x — E^sl. Since the presence 
of a magnetic field at arbitrary direction would overcom- 
plicate the calculations, we shall consider in the following 
the case of a zero magnetic field. The TDGL equation 
analogous to Eq. ((HJ will write in this case: 



(28) 



+ao- (1 - cosgs) ^{k,ky,q,t) = C{k,ky,q,t), 



where the Fourier transformed order parameter 
■>pik,ky,q,t) is given by 



dk f dky dq 



(29) 



.Q-ikx ^-ikyy ^-iqls 



ip{k,ky,q,t) , 



VI. NON-OHMIC CONDUCTION FOR A 
TILTED ELECTRIC FIELD 

In the linear response approximation, the current pro- 
duced by an arbitrarily oriented electric field can be sim- 
ply obtained by superposing the currents generated sep- 



are 



and the noise terms 
as{ak,ky,q,t)C {k',k'y,q',t')) 

" k')5{ky 



delta-correlated such 
= 2T^^kBT{2nf ■ 
5(k - k')S{ky - k'y)5{q- q')5{t - t'). It can be verified 
that with the aid of the Green function in the three 
variables {k,q,t) for Eq. H28|l . the solution for the 
Fourier-transformed order parameter writes: 



'ilj{k,ky,q,t) = Fo / dT({k 



■ exp < -aoFoT 



2eEx 



T, ky, q 



2eE, 



-T, t - t) 



(30) 



2eEr 



^0 f^eEx 



■ exp 



aoTohr eE^sr 
——= — sm — - — cos I qs 
2eE,s h ' ^ 



2 

eE,ST 



Consequently, the general order parameter correlation function that will allow us to calculate the Cooper pair density 
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Temperature [K] 



FIG. 1: (a) Out-of-plane resistivity as a function of temperature in YBa2Cu306+a;, for two values of the magnetic field, at 
several magnitudes of the electric field, when both AL and DOS fluctuation contributions are taken into account. The material 
parameters are given in the text. The arrows show the increasing electric field direction, (b) Same as (a), when only the AL 
part is considered, (c) The out-of-plane AL conductivity, at different magnitudes of the electric field, at B = 5T. (d) The DOS 
contribution to the out-of-plane conductivity, at different electric fields and fixed magnetic field. 



as well as the current density along the x and z axes will write: 



ao Jo 



du\ — 



dk f dk. 



27r / 27r 



dq 

277 



exp 



£'sin6' 



I r. AT^ E sin 2 \ 
•exp ( 2v3foK — — u ] exp 



r sin {j)u cos 0) 
2 p cos 9 



cos (gs — pu cos I 



where Eq and p are defined in Eq. H10() . One should also notice the cut-off condition for the in-plane momentum 
eo{kl + k^)<c. 

We are now able to summarize the results obtained for the Hartree renormalization equation and for the current 
density components as being, respectively, 



ln^+,T rdue-<'-^)-<^y-'lJ'-'^^\ ^d.e-»%f2^/3V^^. 
To Jo V2 Pz J Jo V Eo 



(31) 



o2 



dwwe—iIo^h)^^^^^^^,y (32) 



where Pz = pcos6 = (sa/S/^o) {Ez/Eo), 7a = ^Oc/£.o = V^V"^ is the anisotropy factor, and Iq, Ii and I2 are the 
modified Bessel functions. The current density components can be written in a simpler form if the cut-off procedure 
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is neglected (c oo), namely 



2 7^ /"^ J ^-u(l+r,\-(EJE„)^u^ sin(p^u) {rsui{p^u 



u(e+^)-{EJE„fu^ , sin(p^M) ^ / r sir 



Besides the superconducting fluctuation contribution 
{AjxT Ajz), also the normal state conduction must be 
considered in the total current density, which will be thus 
given by: 



J, = Aj.,[E^,E,) + a'^^E, 
j, = Aj,{E,,E, 



(36) 



if one supposes Ohmic in-plane and out-of-plane normal 
state conductivities cr^^ and cr^, respectively. As it is 
generally the case for an anisotropic conductor, the cur- 
rent is not parallel to the fleld, unless the latter is applied 
along one of the principal axes of the material. Moreover, 
unlike the Ohmic regime, the current densities compo- 
nents Ajx and Aj^ in Eqs. and lO or in Eqs. 
and H35|) depend generally on both E.j. and E^, so that 
the non-Ohmic effect of an arbitrarily oriented electric 
field can not be reduced to the superposition of the non- 
Ohmic currents produced separately by the in-plane and 
out-of-plane field components. 

It must be also mentioned that the superconducting 
fluctuation current ( Aj j;, Aj^) in Eq. (|36|l . as inferred in 
the GL approach, regards only the AL fluctuation pro- 



VII. RESULTS IN THE 3D LIMIT 

A. General case of non-zero magnetic and electric 
fields 

The results obtained in the previous sections for a lay- 
ered superconductor cannot be transformed for the three- 
dimensional case by directly imposing the 3D condition 
s — *■ (or r —^ oo), because in the layered model it was 
assumed that a cut-off in the z direction is not necessary, 
since the out-of-plane momentum q was already conflned 
to the interval [— tt/s, tt/s]. However, in the 3D case, a 
cut-off condition is necessary for all the three momen- 
tum components, {kx,ky,q). Assuming isotropy in the 



xy-plane, the 3D cut-off condition for the total kinetic 
energy will write thus, in the absence of a magnetic field. 



2m 2m 2m 



< c. 
(37) 

If a magnetic field B is applied along the z-axis, the cut- 
off condition will be written though in terms of the Lan- 
dau level rt, as: 



aoh{2n+l) 



2mr 



< cao or h{2n + 1) + g^Coc 



< c. 



(38) 

The kinetic energy ti^q^ /2mc will replace therefore the 
interlayer coupling energy aor (1 — cos gs) /2 from the 
layered case. This can be formally performed by taking 
the cosine function in the small-^/- limit, such as: 



r 

2(' 



COS qs) 



9 ? 

rq^s^ 



1 SOc 



(39) 



which can be interpreted by saying that the 3D behavior 
is approached when the size of the Cooper pairs along 
the z-axis is so large that the peculiarities of the layered 
structure do not play a role any more, meaning that only 
small values of the out-of-plane momentum q are impor- 
tant in the integrations. This regime indeed occurs for 
the layered superconductors in the near vicinity of the 
transition temperature, where the fluctuations acquire 
anisotropic 3D character. 

In the presence of an electric fleld E along the z-axis, 
one must solve now the TDGL equation in the form: 



T^-i d d _ 
° at ^ + ao£n + aoq 



■ipn{k,q,t) = Cn{k,q,t), 
(40) 



so that the solution is foimd to be: 



tpn{k,q,t) ^ Tq I dT(n{k,q-'^^T,t~T) -expl-ToT 



aoEr 



eE 
q —T 



6mc 2toc ^ 

The current density along the z-direction in the 3D case has in the chosen gauge the usual form: 



(41) 
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0.) 



2mc 



(9^ - dz>) {ip {x, y, z, t) ■0* (x, y, z' , t)) 



where 



n>0 

and the correlation function is, according to Eq. H41|l . 



Hk 
2^ 



(42) 



(43) 



{'ip{x,y,z,t)'ip*{x,y,z',t)) = /t 



• exp 



E 



i—{z-z) 

t.Oc 



dq' 

n / 



(44) 



h(2n+l)+q''^<c 

■ exp {~u [e + (2n + 1) ft + g'^ - 2g'p'it] } , 



where 



p ~ ——, — — tj — 



^0 Eq 



(45) 



After careful evaluation of the q' integral and the LL sum, with consideration of the cut-off condition H38|l . one 
obtains eventually the Hartree rcnormalization equation with both the magnetic and electric fields applied along the 
z-direction, as: 



?{3D) 



T , .9 



(3D)y 



du ■ 



2uh e 



-uh 



1 — exp {—2uh) 



c+h sinh \ 2p'u^Jw\l ^ 







2p'u 



,',,2 



, (46) 



3 



(3D) 



(47) 



One should note that kql in Eq. (|47() is the Ginzburg-Landau parameter in the xy-plane (and therefore proportional 
to the Cooper pair mass m in the a;?;-plane, which in turn is proportional to l/^^o); that the g^^'^-' -parameter is in 
fact proportional with the product symmetric with respect to the coherence lengths along the three principal 

axes of the material. 

The general expression of the fluctuation conductivity (the AL contribution) in a 3D anisotropic superconductor, 
when both the electric and magnetic fields are applied along the symmetry axis, will write in turn: 



i-c+h 

Jo 



3/2 



2uhe 



duu^ 7 — „ 

Q 1 — exp {—2uh) 



exp 



4p'2y3 



-ue ■ 



dw u;3/2e-"'" 



cosh 2p'u'^J^ sinh 2p'u'^J^ 



(2p'M2^y^) {2p'u^^. 
The above equation H48(l takes a simpler form if one neglects the cut-off (i.e. for c — ^ oo), namely: 



du 



2uh ( 



32-^71/1^0 Jo 1 ~ 6xp {—2uh) y/u 



exp 



-ue p'^-u^ 

3^ 



r 



(48) 



(49) 



This expression can be also directly inferred from the 
corresponding result of the layered model, if one takes the 
3D limit s — > 0, r oo in Eq. p4|l . In this case p — > 0, 



so that sm{pu)/pu 1, while the sine function can be 
expanded up to the cubic term in u in the argument 
of the modified Bessel function, which becomes thus 
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h {ru/2 - eocE'^u^lQ)-{^r)-^'^ exp (rw/2 -p'^uVs), 
if one takes the asymptotic expression of the modified 
Bessel functions /„(z) ~ e^/\/27rz for large arguments 
z ^ oo. It must be reminded, however, that the general 
expression (|48|l . where the cut-off is considered, can not 
be directly obtained from the corresponding result H13() 
of the layered model, since for the latter the cut-off 
procedure is not applied for the out-of-plane momentum. 



netic field, can be found directly from the correlation 
function (|44|) . by performing the q' integral before the 
LL sum. They write: 



B. Limit cases of vanishing electric or /and 
magnetic field 

The linear response limit for the results 146() and (|48|l , 
i.e. the case of a vanishing electric field and a finite mag- 



arctan - 



I c + e 



A AL|(3D) 



E 



1 

-3/2 



arctan - 



1/2 



(50) 
(51) 



where the LL cut-off index is — {c— h)/2h and e„ = £ + (2n + l)h. 

In the other limit of a vanishing magnetic field {B = 0) but under arbitrarily strong electric field {E > 0), one 
obtains: 



~{3D) 



In 



Tn 



du exp —ue — 



4p'2y3 



dw e 



sinh (2p' 

-uw \ ^ 



(52) 



which is the 3D equivalent of the renormalization equation 117(1 . and 



Act 



AL|(3D) 



SttHQ Jo 



duu^e-''^- 



„'2„3 /"C 



I 



cosh {2p'v?^Jw) sinh (2p'i 
{2p'u'^^f 



{2p'u'^y^y 



, (53) 



which is the 3D equivalent of the fluctuation conduc- 
tivity 1)18(1 from the layered model. Taken in the isotropic 
case (^oc — io), Eqs- and match the correspond- 
ing expressions already found in Ref. 19, where the cal- 
culations were performed in 3D limit in connection to the 
non-Ohmic in-plane conductivity for the layered model. 

If one neglected the cut-off (c — > oo), the r.h.s term in 
Eq. ((52(1 would become divergent, while Eq. 153(1 would 
take the expression 



Act 



zz iNoCut 



320Ffieo' Jo 



du 



-eu—^p u 



already knownMiiSi^ for isotropic bulk supercondutors 
(Coc = Co) and Gaussian fluctuations (i.e. with £ = £). 



Taking further the limit £^ ^ 0, one obtains for the 
Hartree renormalization equation: 



~<3D) 



In 



T 2g(3D)T 



Tn 



' £ arctan 



and the fluctuation conductivity in the z direction. 



(55) 



Act 



AL|(3D) 

l£;=o,s=o 



487rftCo 
3£^/c 



3 arctan 



V£ 
5c3/2 



(56) 



(c + ?)' (c + ?)'_ 

The result ((56() matches thus the expression foundi2L2£ 
for Gaussian fluctuations [e — e) in a 3D isotropic super- 
conductor (Coc = Co)- 
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Equation (|55|l taken for e = gives the relation be- 
tween the mean-field transition temperature Tq and the 
actual critical temperature Tco = 2^cIe=o b=o where the 
superconductivity is attained in the absence of external 
fields, 

To = T,o exp (25(30)^^^^/^^ (57) 

which represents the equivalent of Eq. H21|) from the 
layered model. 



C. Estimation of DOS term 



Supposing that the same proportionality (|24|) between 
the DOS contribution to the fluctuation conductivity 

Acr^jP^ and the Cooper pairs density holds also 

in the 3D case, one can estimate for the different combi- 
nations of electric and magnetic fields. 



A, 



DOS 

zz 



A<T^i''= (E,B)r ' = / du- , „ ,, e-"'^ — / dzwe"""' " ^ , o . (58) 

^ ^' 47rn?o io l-exp(-2u/i) 



(^^1-0 ^ -^^E7--tanW^-l, (59) 



n=0 



A-?°^lS,.^o - (V^- Vfarctan^l) . (61) 



Unlike the layered model, where, for instance, the AL 
contribution H15() was quadratic in the anisotropy param- 
eter r while the DOS one H22|) was linear, the correspond- 
ing 3D fluctuation contributions (|51|l and H59|l are both 
proportional to the ratio ^oc/Co' while the AL one is more 
singular in £. The DOS contribution might thus be of 
less importance in the 3D case with respect to the lay- 
ered model. However, the relations (|58|l . I|59|) . (|60|l and 
(|61|l remain to be confirmed or refuted by a microscopic 
approach. 



D. Arbitrary orientation of the electric field 

The 3D equivalents of Eqs. (|^ and 

(|35|l , valid for an arbitrary orientation of the electric field 
with respect to the layers, can be obtained in a sim- 
pler manner, without having to solve again the TDGL 
equation, by using a special scaling transformation of 
the coordinates and field components that reduces the 
problem to the isotropic case>i2i2S In the isotropic sys- 
tem, the coordinate axes can be in turn freely rotated so 
that the electric field acquires again only one non-zero 
component, for which the solution is already known. For 
the anisotropic 3D model with axial symmetry treated in 
this Section, having the zero-temperature-extrapolated 
coherence lengths ^oc in the ^-direction and ^0 in the xy- 
plane, corresponding to the anisotropy factor 



7a 



?0 



< 1 



(62) 



the scaling transformation of the coordinates and vector 
potentia l^^i™ writes: 



x = x; y = y; z 
This transformation implies the rescaling of the fields as 



z 

7a' 
Az = 7a^z , 



(63) 



Bx = laBx 

Ex = Ex 



By — JaBy] Bz 



B, 



Ey = Ey- Ez=-iaEz- (64) 
and for the order parameter and current density as 



3x — lajxj 



ia\ip\ ; 

3y — lajy\ Jz — 3z 



(65) 



Thus, the anisotropy is removed from the gradient terms 
in the GL equation and reintroduced into the magnetic 
energy term, whose fluctuation can however be usually 
neglected for hard type-II superconductors.'^^ 

In the case of a zero magnetic field and an electric 
field {E3; = EshiO, Ey = 0, i?2 = EcosO) apphed in the 
anisotropic material at an angle 9 with the z-axis, the 
problem can be reduced, according to Eqs. (|64|l . to that 

of an electric field (^E^ = EsinO, Ey — 0^ E^ — EcosO^ 
applied on an isotropic superconductor having the coher- 
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ence length ^o, such as: Taking now into account Eqs. ifS^ . and ifS^ con- 

sidered for ^oc = £.0, together with Eqs. (|^ . one obtains 
E = E (sin^ 61 + 7^ cos^ 9) ^'"^ , eventually 

tans' = —tan 6*. (66) 

7a 



S<3D) 



In- 



To TT 



du exp 



-we — 4 



St^Hoc Jo 



J 2 -"£-4 u-^ 



9 —ue — 4:( ^ \ iL^ 

du u e V ^0 J 



Eo 



dww^/'^e-"'" 



dww^'^e-'''" 



sinh 



dw e 



2\/3- 



S „2 



cosh(^2V3-^u^V^«) sinh (^2a/3;|^m2^ 



(67) 



cosh (^2^3;^ 1*2 Ti^;^ sinh (^2V3;^u27iy 



2\/3-^u2vw 



(§9) 



so that neglecting the cut-off, 



(3D) 



NoCut 



NoCut 



320r?i^oc Jo nA* 



: exp 



; exp 



'£ U - 



-e u — 



[eo^ 
E^j 



(70) 
(71) 



VIII. CONCLUSION 

In this work we have theoretically investigated the non- 
Ohmic effect of an arbitrarily strong electric field on the 
out-of-plane fluctuation magneto-conductivity of a lay- 
ered superconductor. Our framework was provided by 
the Langevin approach to the TDGL equation, and the 
Hartree approximation was used in order to take into 
account the fluctuation interaction. The main general 
results of our work, valid when magnetic and electric 
fields of arbitrary magnitude are applied perpendicular 
to the layers, are the formulae Ullfl for the renormalized 
reduced temperature and, respectively, (fT^ for the AL 
contribution to the out-of-plane fluctuation conductivity, 
as well as the estimation (|26|) for the DOS term. In the 
limit case of a vanishing electric field, the results were 
found to reduce to the expressions already known from 
the linear response approximation fSection IIII|1 . Exten- 
sions of the results have been provided for the case of 
a tilted electric field with respect to the crystallographic 
axes f Section I VI |) . as well as for 3D anisotropic supercon- 
ductors fSection [VII|l . 

In order to illustrate the predictions of the theoretical 
calculations, we have taken as example a typical HTSC 
material, like the optimally doped YBa2Cu306+2:, and 
evidenced the non-linear effect of a strong electric field 



through comparison to the results obtained in the lin- 
ear response approximation, in the presence of a finite 
or of a zero magnetic field. An important fluctuation 
suppression in the AL part of the out-of-plane conduc- 
tivity has been predicted for electric flelds of hundreds of 
V/cm, while the non-Ohmic effect on the DOS contribu- 
tion turned out to be marginal in the same range (Fig. 

m 

So far, the effect of a strong electric fleld on the 
transport properties of HTSC has been investigated 
and proven experimentally only for the in-plane fluc- 
tuation conductivity and only in the absence of mag- 
netic field.2°i2ii22,24 rpj^g difficulties lie in the high dis- 
sipation in cuprates (of the order of GWcm"'^ at elec- 
tric fields of hundreds V/cm) that can increase the sam- 
ple temperature to values where the non-linearity is no 
longer discernable. In this connection, applying short 
current pulses (tens of ns) at high current densities (a 
few MAcm~2)^ in combination with using very thin films 
(under 50 nm thick) in order to enhance the heat trans- 
port to the substrate)^ seems to be a better alternative 
to the dc and ac measurements. This procedure, if ap- 
plied on vicinal films in order to access also the c-axis 
conduction, will probably allow for the necessary accu- 
racy necessary to detect the non-Ohmic behavior also for 
the out-of-plane fluctuation conductivity. 
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